Effect of hydrogen adsorption on the quasiparticle spectra of graphene 
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We use the non-interacting tight-binding model to study the effect of isolated hydrogen adsorbates 
on the quasiparticle spectra of single-layer graphene. Using the Green's function approach, we 
obtain analytic expressions for the local density of states and the spectral function of hydrogen- 
doped graphene, which are also numerically evaluated and plotted. Our results are relevant for 
the interpretation of scanning tunneling microscopy and angle-resolved photoemission spectroscopy 
data of functionalized graphene. 

PACS numbers: 73.20.Hb, 73.22.Pr, 78.67.Wj 
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The adsorption of hydrogen on graphene is of funda- 
mental as well as applied significance. Current inter- 
est in hydrogenated graphene is due to the possibility 
of opening a band gap in the semimetallic graphene, 
which is essential for applications in electronics. This 
has recently prompted several research groups to experi- 
ment with hydrogenated graphene and to study its elec- 
tronic and transport properties Calculations based on 
density-functional theory (DFT) gave the motivation to 
search for hydrogenated derivatives of graphene, since 
they showed that graphane, i.e., the fully hydrogenated 
compound, is both a stable 2D hydrocarbon and a band 
insulator.— Earlier DFT calculations had studied the ef- 
fect of hydrogen adsorption on the electronic structure 
of graphene, using a 4 x 4 regular arrangement, i.e., an 
H:C ratio of 1:32, and found a band gap opening and 
the appearance of a spin-polarized dispersionless midgap 
band.— 

In this work, we analyze the effect of a single hydro- 
gen impurity on the electronic structure of graphene. As 
we show below, the analysis can be extended also to cal- 
culate some properties of randomly distributed H adsor- 
bates at a low-density coverage. A hydrogen atom ad- 
sorbed on graphene forms a covalent bond atop a carbon 
atom, and changes the orbital state of the host atom 
from sp 2 to sp 3 . This is expected to have a character- 
istic effect on the local density of states (LDOS), and 
to produce Friedel oscillations in the surrounding charge 
density. These properties can be mapped out by scan- 
ning tunneling microscopy (STM) and, indeed, hydrogen 
monomers and dimers on graphene have been imaged by 
this techniquei£i£ As the graphene sample is exposed to 
a hydrogen plasma, we initially expect the adsorption of 
a low concentration of randomly distributed H atoms, 
and in later stages more complicated and perhaps regu- 
lar adsorption arrangements. We can enumerate several 
changes that can be expected in the initial stages of ad- 
sorption to occur in the quasiparticle spectra. First, due 
to scattering from the H impurities there is an effect of 
line broadening in the momentum-resolved spectral func- 
tion. Second, H adsorption breaks the sublattice symme- 
try which can lead to a band gap opening at the Dirac 



point. Third, since hydrogen impurity on graphene acts 
as a resonant scatterer, a midgap resonance appears in 
the density of states^ The spectral function is accessible 
by angle-resolved photoemission spectroscopy (ARPES), 
and the above changes have already been obeserved in 
experiments on hydrogenated graphene^— 

Theoretical studies of the electronic and transport 
properties of hydrogenated graphene have been based 
on a tight-binding model which describes the ir bands 
of graphene and an additional orbital for the H 
adsorbateJ^— The purpose of this Brief Report is to cal- 
culate the local density of states and the spectral function 
of sparsely hydrogenated graphene within the framework 
of this tight-binding model. As in earlier theoretical stud- 
ies of the effect of impurities on the spectral properties 
of graphene, a Green's function analysis can be used to 
obtain the desired quantities.— >±2r— However, in contrast 
with most previous works which have considered substi- 
tional impurities, we deal with adsorbate impurities. 

We begin our derivation by defining the non- 
interacting tight-binding model that describes the 7r 
bands of graphene with an additional adsorbed hydro- 
gen atom a a 10 i 12 



rl — Ho + tL\ + Hi, 
where Ho describes pure graphene, 
Ho = -t^c\c 3 , 

(id) 

Hi describes an isolated hydrogen atom, 
Hi = Edd^d, 



(1) 



(2) 



(3) 



and Hi describes the hybridization between graphene 
and hydrogen adsorbed on an arbitrary host site spec- 
ified by index I, 



Hi 



V(c]d + 



(4) 



Interactions involving spin have not been included in this 
model, and the spin degree of freedom has been omit- 
ted for simplicity. The energy parameters which have 
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FIG. 1. (Color online) Graphical solution of e — £d — 
V 2 Rego(e) = showing the positions of the resonance just 
to the left of Dirac point, and two localized states just below 
and above energy bands, respectively. 



been obtained from fits to first-principles band struc- 
ture calculations are given by t = 2.6 eV, V = —2t and 
e d = -t/16ii 

The local density of states at a given site m is given 

by 



p(m;s) = lmG(m, m; e), 

7T 



(5) 



where G(m, m; e) is a matrix element of the Green's func- 
tion g(e) =(e-H + i0 + ) _1 , while the total DOS of the 
system is given by 



p(e) = ImTr£(e). 

7T 



(6) 



Therefore, we must calculate the Green's function cor- 
responding to Eq. (JTJ. The crucial step is to consider 
the graphene or the hydrogen atom as a separate entity, 
and incorporate the presence of the other system in its 
Hamiltonian in terms of a self-energy. 

The self-energy in a device, is generally given by VcjrV^ 
where V is the coupling Hamiltonian and gn is the subset 
of the Green's function of the reservoir which makes the 
interface with the device^ Thus instead of Eq. ©, we 
can describe a hydgrogen atom in contact with graphene 
by 



Hi = [e d + V 2 g (e)] dU, 



(7) 



where the additional term is the self-energy induced by 
graphene and go(e) is the diagonal matrix element of the 
Green's function of graphene which describes the host 
site. On the other hand, for graphene in contact with 
hydrogen, the Hamiltonian becomes^ 



n = -tJ24 c i 

(i,j) 



V 2 



Ed 



(8) 



where the additional term is the self-energy induced by 
hydrogen at the host site, which has the general form 
if we note that 1/(e — Ed) is the Green's function of an 
isolated hydrogen. We can now calculate the LDOS from 
the Green's functions corresponding to Eqs. ([7]) and (JSJ, 
respectively. From Eq. it follows that the LDOS at 
the hydrogen site is given by 



Pn(s) 



- — Im — 

7T £ - 



1 



£d - V 2 g (e)' 



(9) 



Equation jSJ) is now similar to the Hamiltonian of a 
substitutional impurity, with an energy- dependent on- 
site parameter, so the corresponding Green's function can 
be obtained exactly frorn^i 



Go = Go + GoTGo 
where the i-matrix is given by-iS 

r. 



- £d- V 2 g a {e) 



c\ci. 



(10) 



(11) 



We can find the positions of the resonances and localized 
states by setting the denominators in Eqs. © and (fTTj) 
equal to zero, 



e-e d ~V 2 Reg (£) = 0. 



(12) 



A graphical solution of Eq. (TT21 is shown in Fig. Q] The 
position of a solution just to the left of the Dirac point 
implies a sharp resonance since Im go , which measures 
the width of the resonance, is vanishingly small. It is to 
be noted that the strong hybridization between H and 
C, parameterized as V, is responsible for the resonance 
being much closer to the Dirac point than the on-site 
energy, Ed- The other solutions of Eq. (fl2"|) correspond to 
two localized states with energies just below and above 
the band limits, respectively. It is interesting to compare 
this picture with the case of a vacancy which is also a 
resonant scatterer. For a vacancy, the equation to solve 
would change to 1/U = Remote) where U would now be 
a large positive constant representing the on-site energy 
of a vacancy. There will again be a resonance state near 
the neutrality point, but only one localized state outside 
the bands at a large positive energy, (cf. Ref. Ibfl ) 

From Eq. (fTU)) . it is easy to show that the LDOS at a 
given site of graphene is given by 



/ \ 1 T (e - £d).9o + V 2 (g m ig h 

p{m; e) = Im — 

7r e- e d - V go 



9l) 



(13) 



where g m \ = (m|C?o^), and for the honeycomb lattice we 
can use g m i = gi m . In particular, for the host site m = I, 
and with gu = go, we find 



/ t \ i T (E-£d)go 

p(l- s) = Im — r -. 



(14) 



When the site m is very far from the host site I, g m i — > 0, 
and p(rn;e) — > (— 1 / 7r)hngo(e) = Pq(e), i.e., it approaches 
the LDOS of pure graphene. 



3 



20 




FIG. 2. (Color online) Local density of states of graphene with a single hydrogen adsorbate. (Left panel) The subplots show 
LDOS at the H site, the host site, Co, and its four nearest neighbors, denoted by C1-C4, as well as the change in total DOS 
of graphene, and the DOS per carbon atom for 0.05% hydrogen coverage. (Right panel) LDOS at the H site (at center) and 
nearby graphene sites of a rectangular sample of 6 x 6 unit cells, for an energy near the resonance ss — 0.03eV. 



Evaluation of the LDOS at any lattice site by Eqs. (fT3")) 
and ([TJ]) requires numerical values of the lattice Green's 
functions go and gi m , which can be conveniently calcu- 
lated by the expressions given by Horiguchi in terms of 
elliptic functions^ In Fig. [5J we have plotted the LDOS 
at the host site and its four nearest neighbors, which show 
that the resonance is small on the sites of the sublattice of 
the host, and large on the other sublattice, which seems 
to be a common property that is shared by vacancies 
also Figure [5] shows that at sites farther from the im- 
purity the LDOS approaches that of pure graphene in an 
oscillatory fashion. 

From Eq. (fTU)) it follows that the change in the total 
DOS of graphene is given by 



Ap- 



1. 



-Im 



dgo 



V 



de e - E d - V 2 g 



(15) 



To calculate the DOS corresponding to Eq. (JXJ) , since we 
need a finite number of impurities to make the effect fi- 
nite, we assume to have N carbon atoms in the graphene, 
and Ni <^ N well isolated hydrogen adsorbates. We can 
then write the DOS per carbon site as 



p(e) = p (e) + n t [Ap(e) + pn(s)], 



(16) 



where po is the DOS of pure graphene, and n$ = Nt/N 
is the ratio of hydrogen to carbon atoms. A plot of pn 
is shown in Fig. [2] where the expected resonance features 
can be observed. We have also plotted the change in the 
total DOS, Ap, and the DOS for graphene with a tiny 
amount of 0.05% of hydrogen. We note that the total 
DOS is in agreement with those obtained by numerical 
calculations^ It must be remarked that the values of 
Ap can become negative while the total DOS must al- 
ways remain positive. The explanation why this is not 
unphysical is that the measurable quantity, which is the 
total DOS, remains positive, because it is of O(N) while 



the change in it due to the presence of a single H adsor- 
bate is of O(l). Also shown in Fig. rj is the LDOS at 
the sites of a rectangular sample of 6 x 6 unit cells of 
graphene, at an energy near the resonance. This result 
shows the threefold symmetry and fluctuations of LDOS, 
and can be related to STM images. 

It is easy to extend the above results to find the effect 
of a low density of H adsorbates on the spectral function, 
given by 



A(k,E) = 



-2ImS(k,e) 



[e - £ k - Re E(k, e)] 2 + [Im £(k, e)] 2 ' 



(17) 



where S(k, e) is the self-energy due to the presence of 
hydrogen impurities and £k is the energy dispersion of 
graphene, 



£k = ±t 



2 cos K-a + 2 cos 



— — cos — =— — 

2 2 



1 1/2 



(18) 

where a = 2.46 A is graphene lattice constant. For low 
concentrations of H adsorbates, the average t-matrix ap- 
proximation (ATA) applies and, noting that here we have 
no explicit k dependence, we can write 



E(e) 



n,V 2 



e- £ d - V 2 g (£) ' 



(19) 



The energy and momentum widths of the spectral func- 
tion at a given energy near the Dirac point can be ap- 
proximated as Aek = ft/Tk = — 2Im£(ek); and KAk — 
— (2/vp ) Im £(sk)> respectively, where r k is the scatter- 
ing time and vp is the Fermi velocity. Energy and mo- 
mentum widths, which are proportional to rij at low con- 
centrations according to Eq. (|19p . can be deduced from 
ARPES energy distribution curves (EDC) and momen- 
tum distribution curves (MDC), respectively, and can be 
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FIG. 3. (Color online) Spectral function of graphene along the 
TKM path of its Brillouin zone. (Left panel) Clean graphene. 
(Right panel) Graphene with 2% hydrogen adsorption. Gen- 
eral line broadening, appearance of a localized state just be- 
low the band, and changes similar to a gap opening near zero 
energy are the noteworthy features. 

used to estimate the density of randomly distributed H 
adsorbates using these results. In Fig. [31 we plot the spec- 
tral function for pure graphene and for graphene with a 
relatively large 2% impurity concentration to make the 



effects large enough to be more clearly visible. Ourplots 
compare well with the experimental results of Ref. [j. We 
can observe the effect of spectral broadening and the fin- 
gerprints of the localized state below the band boundary 
at —7.8 eV. Less trivial are the changes observed around 
the Dirac point, i.e., the zero energy at K point, where 
both sublattice symmetry breaking and the presence of 
the midgap resonance can play a role. Even within the 
simplest approximation that we have used to relate the 
many-impurity to the single-impurity scattering, we see 
a quasi band-gap opening and the appearance of new 
quasilocalized states within the gap near zero energy^^ 

In conclusion, we presented analytic expressions for the 
LDOS of graphene with a single hydrogen atom adsor- 
bate, as well as for the spectral function of graphene with 
a low density of randomly distributed hydrogen adsor- 
bates. Our plots of the spectral function explain the fea- 
tures of quasi band gap and the spectral line broadening 
observed in angle-resolved photoemission spectra of hy- 
drogenated graphene. Our expressions can be used to 
estimate the concentration of hydrogen adsorbates from 
measured spectral linewidths. 
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